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Abstract 

Analytical models were developed to model the 
heat transfer through high-temperature fibrous 
insulation used in metallic thermal protection systems 
on reusable launch vehicles. The optically thick 
approximation was used to simulate radiation heat 
transfer through the insulation. Different models for 
gaseous conduction and solid conduction in the fibers, 
and for combining the various modes of heat transfer 
into a local, volume-averaged, thermal conductivity 
were considered. The governing heat transfer equations 
were solved numerically, and effective thermal 
conductivities were calculated from the steady-state 
results. An experimental apparatus was developed to 
measure the apparent thermal conductivity of insulation 
subjected to pressures, temperatures and temperature 
gradients representative of re-entry conditions for 
launch vehicles. The apparent thermal conductivity of 
an alumina fiber insulation was measured at nominal 
densities of 24, 48 and 96 kg/m 3 . Data were obtained at 
environmental pressures from 10 4 to 760 torr. with the 
insulation cold side maintained at room temperature 
and its hot side temperature varying up to 1000°C. The 
experimental results were used to evaluate the 
analytical models. The best analytical model resulted 
in effective thermal conductivity predictions that were 
within 8% of experimental results. 


Nomenclature 

A ratio of parallel to total heat conduction 
c specific heat 
D, fiber diameter 
d gas collision diameter 
e specific extinction coefficient 


f solidity ratio 
G constant 
k thermal conductivity 
K b Boltzmann constant 
L thickness 
L c characteristic length 

m number of control volumes for the numerical 
model 

P pressure 
Pr Prandtl number 
q" heat flux 
T temperature 
t time 

x spatial coordinate 

Greek Symbols 
a accomodation coefficient 
(3 extinction coefficient 
AT temperature difference across test specimen 
£ emittance 
y specific heat ratio 
X mean free path 
p density 

a Stefan-Boltzmann constant 
T 0 optical thickness 

Subscripts 
a apparent 
c conduction 
e effective 
g gas 
r radiation 
s solid 
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Introduction 


Metallic thermal protection systems are being 
considered for insulating major portions of a new 
generation of reusable launch vehicles (RLV). The 
main function of the thermal protection system (TPS) is 
to maintain the vehicle structural temperature within 
acceptable limits during re-entry aerodynamic heating. 
The metallic TPS consists of a metallic shell panel 
fabricated from high-temperature metallic alloys and is 
mechanically attached to the substructure. The shell is 
filled with lightweight, non-load-bearing, high- 
temperature fibrous insulation. 1,2 Insulation for current 
metallic TPS concepts experiences environmental 
pressure from 10 2 to 760 torr, while the hot surface of 
the insulation is exposed to temperatures as high as 
1000°C. 

Heat transfer through a fibrous insulation involves 
combined modes of heat transfer: solid conduction 
through fibers, gas conduction and natural convection 
in the space between fibers, and radiation interchange 
through participating media. The relative contributions 
of the different heat transfer modes vary during re- 
entry. Radiation becomes more dominant at high 
temperatures and with large temperature gradients 
through the insulation, while gas conduction and 
natural convection contributions are minimal at low 
pressures and become more significant with increasing 
pressure. The complex coupling of the heat transfer 
modes makes the analysis and design of high- 
temperature insulation difficult. 

Heat transfer through fibrous insulation has been 
investigated by various researchers in the last 40 years, 
but the majority of the work has been limited to small 
pressure and temperature ranges. Verschoor, et al.,' 
experimentally and theoretically studied heat transfer 
through fibrous insulation over a wide pressure range, 
but only for temperatures up to 150°C. Larkin and 
Churchill 4 performed an experimental and theoretical 
study of heat transfer in fibrous insulation at 760 torr 
and temperatures up to 400°C. Hager and Steere 5 
modeled the heat transfer through fibrous insulation and 
compared it with experimental results at room 
temperature at a pressure of 10 4 torr. Bankvall 6 
provided a theoretical simulation of heat transfer 
through building insulation and compared it with 
experimental results obtained at pressures from 10 to 
760 torr, and at temperatures up to 80°C. Pawel, et al., 7 
modeled the heat transfer through fibrous insulation and 
compared it with experimental results for pressures 
between 40 and 760 torr and temperatures up to 930°C. 
Williams and Curry 1 have provided a theoretical 
formulation of heat transfer through rigid fibrous 
insulation. Tong, et al., 9, 10 modeled heat transfer 
through insulation at atmospheric pressure with 
temperatures up to 150°C. 


The goal of the present research is to evaluate 
various analytical models for calculating effective 
thermal conductivity in fibrous insulation by comparing 
predicted values with measured values of apparent 
thermal conductivity. To facilitate this evaluation, a 
thermal conductivity apparatus was developed for 
measuring apparent thermal conductivity of high- 
temperature fibrous insulation at pressures from 10 4 to 
760 torr. and temperatures up to 1000°C. This paper 
describes the analytical formulation for heat transfer in 
fibrous insulation, the experimental apparatus for 
measuring the apparent thermal conductivity, and the 
comparison of measured and predicted thermal 
conductivities. 


Analytical Model 

The heat transfer mechanisms in fibrous insulation 
include solid conduction, gas conduction, natural 
convection, and radiation. Natural convection is 
significant at environmental pressures above one torr. 
In the present study, the insulation is assumed to be 
confined between two large horizontal plates, with the 
bottom plate having a lower surface temperature than 
the top plate. This assumption eliminates natural 
convection heat transfer. 

A continuum model for conservation of energy in 
one-dimensional heat transfer in an insulation by 
conduction and radiation yields the partial differential 
equation: 1 1 
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d_ 

dx 
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Two boundary conditions and an initial condition are 
needed to complete the problem statement. The 
insulation is assumed to be at a uniform initial 
temperature, and constant temperature boundary 
conditions are assumed to simulate the hot and cold 
surfaces of the steady-state thermal conductivity 
apparatus. 


Radiation 


In its most general form, the radiant heat flux, q", 
is defined by a system of simultaneous differential and 
integral equations. Various techniques have been used 
for addressing the radiation component of heat transfer 
through insulation. The optically thick approximation 
is valid for optical thickness much larger than unity. 

The optical thickness is defined as the ratio of the 
characteristic length to photon mean free path. Because 
the optical thickness of the insulation samples used in 
this study varied from 13 to 86, the optically thick 
approximation was used. By assuming a gray medium 
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and utilizing the optically thick approximation, the 
radiant flux can be reduced to:" 


q, =- 


16a T 3 3T 
3f) 5x 


( 2 ) 


This formulation is valid in the interior of the medium 
but not at the boundaries, because it does not include 
any terms for radiation from the bounding surfaces. 
Despite this limitation, it is assumed that the 
formulation is applicable over the entire domain. A 
radiant thermal conductivity can be defined as: 


k =i^T> 
' 3|J 


(3) 


The extinction coefficient, (1. is the inverse of the mean 
penetration distance of radiation in absorbing and 
scattering medium, and is given by: 

|3 = ep (4) 

where e is the specific extinction coefficient. The 
variation of specific extinction coefficient of alumina 
fiber with temperature given by Keller and 
Blu men berg' was used in the present study. 


Gas Conduction 


Another important heat transfer mechanism in a 
fibrous insulation is the gas conduction in the space 
between fibers. The variation of thermal conductivity 
with temperature for various gases is well established. 
Two models used for the variation of gas thermal 
conductivity in porous media with environmental 
pressure are discussed. Using the temperature jump 
theory for heat transfer in rarefied gas flow, the gas 
thermal conductivity in porous media is given by: 13 
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with k g being the temperature-dependent gas thermal 

conductivity at atmospheric pressure. Marcussen' 4 has 
used the same formulation for modeling variation of gas 
thermal conductivity in fibrous insulation. The mean 
free path, the mean distance traveled by molecules 
between collisions, is given by: s 


1 = 


k b t 

V27td;P 


The characteristic length, L., is defined as: 36 


(5b) 


L =2 A. (5c) 

c 4 f 

The solidity ratio, f, is defined as the ratio of the density 
of fibrous insulation to the density of fiber parent 
material. As pressure approaches 760 torr, the ratio of 
mean free path to characteristic length in Eq. 5a 
approaches zero, resulting in k g = k . The 
accomodation coefficient, a, is a measure of the 


average efficiency of the energy exchange between gas 
molecules and fibers. Various values for the 
accomodation coefficient between zero and unity were 
utilized in this study to determine the optimal value. 


Various researchers 3 6 15 have used the following 
formulation for the variation of gas thermal 
conductivity with pressure for fibrous insulation: 



( 6 ) 


This formulation was evaluated in this study and was 
found to be unsatisfactory. The numerical predictions 
obtained using this equation rose faster with increasing 
pressure compared to both experimental results and 
numerical results obtained using Eq. 5a. 


Solid Conduction 


The solid conduction through fibers has been 
modeled differently by various researchers. The model 
used in this study was: 5 

k=fV (7) 

with k being the temperature-dependent thermal 
conductivity of the fiber parent material. Some 
researchers have used: 15 

k„ = GkjD” (8) 

where the constant G is used for dimensional 
consistency. Using a value of unity for the constant G, 
Eq. 8 was evaluated in the present study and found to 
provide comparable results with Eq. 7 with the value of 
the exponent, n, set equal to 0.5. The solid conduction 
model should implicitly depend on the solidity ratio, f, 
so that k = k when f = 1, and k =0 when f=0. Based 
on this physical reasoning it was decided to use Eq. 7 in 
this study. 

Combined Thermal Conductivity 


Heat transfer through a porous medium involves a 
complex combination of heat transport mechanisms 
which varies with temperature and pressure. For the 
analysis of the macroscopic heat flow through 
heterogeneous media, various empirical and statistical 
approaches have been used to model local volume- 
averaged properties. 16 The two models used in this 
study were: 

k = f(f 2 k’) + (l-f)k g +k r (9) 


k = A[fk;+(l-f)k g ] + 

(1- A) — ^ +k 

(l-fjkj+fk, r 


( 10 ) 
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Eq. 9 models solid and gas conduction in parallel, 17 
while Eq. 10 assumes that solid and gas conduction are 
both in parallel and series, with A being the ratio of 
heat conduction that occurs in parallel. 7 Setting A=f 2 
makes the solid conduction term consistent with Eq. 7. 
At very low pressures where gas conduction is 
negligible both models then yield k = f V + k r . 

Solution Procedure 


By using the optically thick approximation to 
model radiation, the conservation of energy equation is 
simplified to: 


3T 


PC 3t dx 


f dT) 

k— (11) 

V dx J 

The local volume-averaged thermal conductivity, k, 

(Eq. 9 or 10) was incorporated in a 1-D finite volume 
numerical model of the governing heat transfer 
equation given in Eq. 1 1 . The Crank-Nicolson implicit 
algorithm was used to obtain solutions to the transient 
diffusion equation. The constant temperature boundary 
conditions and the environmental pressures used in the 
program corresponded to the test conditions. 


The numerical solution was marched in time until 
steady state conditions were achieved. At steady-state 
conditions, the effective thermal conductivity was 
calculated from the individual thermal conductivities 
for each control volume, by assuming that control 
volumes represent thermal resistances in series: 


k = m 


f m 1 ^ 

if 

i=i k 


( 12 ) 


Thermal Conductivity Apparatus 

An experimental apparatus was developed to 
measure the apparent thermal conductivity of insulation 
under steady-state conditions for comparison with 
analysis. The apparatus was based on the ASTM C- 
201 standard. The main difference between this 
apparatus and the ASTM standard is that thin foil heat 
flux gages are used for measuring the heat flux in the 
present apparatus, while three guarded water 
calorimeters are used in the ASTM standard. The 
apparatus provides apparent thermal conductivity at 
environmental pressures from 10 4 to 760 torr, with 
specimen cold side temperature maintained around 
room temperature, and specimen hot side temperatures 
up to 1000°C. 

A schematic of the apparatus is shown in Figure 1 . 
The apparatus consists of a water-cooled plate and a 
ceramic board picture frame shown in Figure 2. The 
water-cooled plate is made of aluminum and is 304.8 x 


304.8 x 25.4 mm, and has seven coolant passages 
across its length. A picture frame made of refractory 
fiber ceramic board is set on the water-cooled plate, 
with the insulation sample to be tested placed inside the 
picture frame. The ceramic board picture frame is 50.8 
mm wide and its outside dimensions are 304.8 x 304.8 
mm, resulting in a test sample size of 203.2 x 203.2 
mm. For the tests reported here, the ceramic board 
picture frame, and thus the insulation sample, was 13.4 
mm thick. 

The apparatus uses a ceramic-fiber radiant heater 
which is 304.8 x 304.8 mm, and consists of rows of 
iron-chrome-aluminum heater wires. The heater 
surface can be heated to 1 100°C. Since the heater does 
not provide uniform spatial heating, it is not used to 
heat the test samples directly. Rather, a septum plate 
made of Inconel™ 625 with dimensions of 304.8 x 
304.8 x 6.35 mm is used to provide a uniform 
temperature boundary for the test specimen. The 
septum plate is placed directly on top of the test sample. 

The overall assembly of the apparatus includes 
setting the refractory fiber ceramic board picture frame 
on the water-cooled plate, installing the insulation 
sample on the water-cooled plate, and laying the 
septum plate on the sample and picture frame. The 
radiant heater is suspended from the carbon steel frame 
approximately 38 mm above the septum plate. The 
whole assembly is insulated with 25.4-mm-thick 
refractory fiber ceramic board, with additional 
insulation set on top of the heater, to minimize heat 
losses. The apparatus is placed in a 1.5-m-diameter, 
1.5-m-long vacuum chamber capable of providing test 
pressures from 1 x 10 4 to 760 torr. All tests were 
conducted in a gaseous nitrogen environment. Gaseous 
nitrogen was selected instead of air to prevent 
introduction of water vapor into the chamber, thus 
significantly reducing the time required for pumping 
down to very low pressures. 

Instrumentation 

The water-cooled plate is instrumented with nine 
thin film heat flux gages. Each heat flux gage is 25.4- 
mm long and 19.05-mm wide, with a nominal thickness 
of 0. 1524 mm. The gages are thermopiles encapsulated 
in polyimide film, producing a voltage directly 
proportional to the impinging heat flux with a nominal 
sensitivity of 1.05 x I O' 6 V/[W/m 2 ]. Each heat flux gage 
also utilizes a type T (copper/constantan) thermocouple 
for surface temperature measurement. The heat flux 
gages are bonded to the top of the water-cooled plate 
with their lead wires buried underneath the plate top 
surface. Two epoxy filled trenches are made in the 
plate for subsurface burying of the heat flux gage lead 
wires. 
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A polyimide sheet 0. 127-mm thick covers the plate 
at locations not covered by heat flux gages. The 
polyimide sheet is bonded to the plate to provide a 
uniform surface on the water-cooled plate. The top 
surface of the polyimide sheet is spray painted using a 
flat black paint with an emittance of 0.92 throughout 
the infrared spectrum. 

The septum plate is instrumented with 23 metal 
sheathed type K (nickel-chromium/nickel-aluminum) 
thermocouples. The thermocouples were installed on 
the top surface of the septum plate (opposite side from 
the insulation sample) so that uniform contact exists 
between the septum plate and insulation test sample. 
The septum plate was oxidized in an oven at 1000°C 
for 6 hours after the installation of the thermocouples. 
The emittance of oxidized inconel has been reported to 
be 0.85. 19 

The thermocouple and heat flux gage data are 
collected using a personal computer based data 
acquisition system. The thermocouple data are 
converted to temperature using look-up tables by the 
data acquisition software. Raw voltages from the heat 
flux gages are converted to heat flux by applying a 
manufacturer-suggested temperature correction to the 
raw data and then using the manufacturer' s linear 
calibration for heat flux versus voltage for each sensor. 

Data Analysis 

The apparent thermal conductivity of the sample is 
calculated from Fourier’s law: 


Only the data from the central 127 x 127 mm section of 
the test set-up, referred to as the metered region, are 
used for calculating the apparent thermal conductivity. 
There are five heat flux gages on the water-cooled 
plate, and 15 thermocouples on the septum plate in the 
metered region. The apparent thermal conductivity of 
the sample at each of the five heat flux gage locations 
in the metered region is calculated. Then, the data are 
averaged to obtain the average apparent thermal 
conductivity. 


Uncertainty Analysis 

An uncertainty analysis was performed to obtain 
error estimates for the experimentally determined 
apparent thermal conductivity. The procedure for 
calculating the apparent thermal conductivity bias and 
precision uncertainties is that specified by Coleman and 


Steele. 211 The overall uncertainty in the estimation of the 
apparent thermal conductivity is obtained by combining 
the precision and bias errors using the root-sum-square 
method. 

The bias error for the thermocouples was 
determined by performing a comparison calibration in a 
controlled temperature oven/bath with a National 
Institute of Standard and Technology (NIST) traceable 
thermometer. The bias error for heat flux sensors was 
determined by calibration against a high intensity total 
irradiance standard (solar constant lamp). The precision 
error for the thermocouples and heat flux gages was 
determined from the experimental measurements, by 
calculating the standard deviations of the temporal 
variations of each measured quantity. 

The bias, precision, and overall uncertainties for 
apparent thermal conductivity are presented in Table 1. 
The uncertainties are presented as a percentage of 
uncertainty with respect to calculated thermal 
conductivity. The precision uncertainties were obtained 
from measurements on the insulation sample at 24 
kg/m 3 , and corresponded to measurements at pressures 
of 1x10 4 , 1x10 2 , and 760 torr, and nominal average 
sample temperatures of 50, 300, and 500°C. The 
overall uncertainty varied from 5.5 to 9.7%. The 
highest uncertainties were at the lowest sample average 
temperature. 


Table 1. Overall uncertainty estimates for effective 
thermal conductivity. 


Pressure 

(torr) 

Average 

Temperature 

(°C) 

Bias 

Error 

(%) 

Precision 
Error (%) 

Overall 

Error 

(%) 

0.0001 

51.5 

5.86 

7.68 

9.66 

0.01 

51.5 

5.86 

7.29 

9.35 

760 

51.1 

5.96 

4.73 

7.61 

0.0001 

303.3 

3.98 

3.87 

5.55 

0.0105 

303.4 

3.98 

3.87 

5.55 

755 

302.8 

3.97 

3.85 

5.53 

0.0001 

505.0 

3.91 

3.82 

5.47 

0.0114 

505.2 

3.91 

3.84 

5.48 

474 

506.2 

3.91 

3.82 

5.47 


To provide an independent assessment of 
uncertainty, the apparent thermal conductivity of a 
fumed silica board. Standard Reference Material (SRM) 
1459 from NIST was measured using the current 
apparatus. The thermal conductivity of this SRM has 
been reported by NIST to be 0.021 W/mK at 24°C. The 
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apparent thermal conductivity of the SRM sample was 
also measured using the guarded hot plate technique 21 at 
Holometrix, Inc. at average temperatures up to 300°C. 
The results of the measurements are provided in Table 
2. The specimen average temperature, the measured 
apparent thermal conductivity, and the corresponding 
thermal conductivity from the guarded hot plate data 
are presented. The percentage difference between 
measurements using the current apparatus and the 
guarded hot plate data are also presented. It can be seen 
that the data from the current apparatus are within 5.5% 
of the guarded hot plate data. 


Table 2. Comparison of effective thermal conductivity 
measurements on NIST Standard Reference Material 
1459 using the present apparatus and the guarded hot 
plate technique 


Average 

Temperature 

(°C) 

Thermal 

Conductivity 

[present 

apparatus] 

(W/mK) 

Thermal 
Conductivity 
[guarded hot 
plate] 
(W/mK) 

Percent 

difference 

(%) 

53.4 

0.022 

0.0215 

2.3 

53.9 

0.022 

0.0215 

2.3 

95.4 

0.021 

0.0222 

5.5 


Results and Discussion 

The numerically predicted thermal conductivities 
are discussed. The apparent thermal conductivity of of 
Saffil™, an alumina based fibrous insulation was 
measured at densities of 24, 48, and 96 kg/m 3 . Tests 
were conducted with the septum plate temperature set at 
nominal temperatures of 95, 260, 425, 590, 760, 870, 
and 980°C, and the environmental pressure controlled 
to nominal values of lxlO 4 , I x I O ', 1x10 2 , 0.1, 0.5, 1, 5, 
10, 100, and 750 torr. The water-cooled plate was 
maintained at room temperature. All tests were 
conducted in a gaseous nitrogen environment. The 
measured thermal conductivities were compared with 
predicted effective thermal conductivities. 

Numerical Results 

The analytical methods described previously were 
used to numerically calculate effective thermal 
conductivities. The radiation model given in Eq. 3, the 
gas conduction model given in Eq. 5a and the solid 
conduction model given in Eq. 7 were used in the 
analysis. The two models for combining the various 
heat transfer mechanisms into local volume-averaged 
thermal conductivities given in Eq.'s 9 and 10 were also 


used. The effective thermal conductivity was 
calculated from local thermal conductivity using Eq. 

12 . 

The various heat transfer mechanisms are coupled, 
and the net contributions of each component to the 
effective thermal conductivity can not be independently 
assessed. The contribution of each mode transfer mode 
to the local volume-averaged thermal conductivity 
calculated from Eq. 9 for each control volume in the 
numerical model could be assessed from steady state 
numerical results. The results obtained with the hot 
side and cold side temperatures at 1000°C and 20°C, 
respectively, and environmental pressures of 10 4 and 
100 torr are presented in Figure 3. Data presented are 
for an insulation sample with a thickness of 13.4 mm 
and a density of 48 kg/m 3 . The solid conduction, gas 
conduction and radiation contributions to the local 
thermal conductivity and the local thermal conductivity 
are presented as a function of non-dimensional through- 
thickness location in insulation. The non-dimensional 
through-thickness locations of 0 and 1 correspond to 
boundary temperatures of 20°C and 1000°C. Solid 
conduction is almost negligible. Gas conduction is 
negligible at 10 4 torr. At 100 torr, gas conduction has a 
larger contribution to thermal conductivity compared to 
radiation for through-thickness locations up to 0.25. 
This trend reverses at through-thickness locations larger 
than 0.25, with the radiation becoming more dominant 
with increasing temperature. 

Comparison of Numerical and Experimental Results 

The numerical and experimental data are first 
compared at low pressures, where gas conduction is 
negligible. This enables evaluation of models for solid 
conduction and radiation. Then, the numerical and 
experimental results are compared at higher pressures 
where gas conduction is significant. 

Thermal Conductivity in Absence of Gas 
Conduction 

The analytical predictions were compared to 
experimental results obtained at P = 10 4 torr. At this 
pressure solid conduction and radiation were the only 
heat transfer mechanisms. The local volume-averaged 
thermal conductivity using both Eq.’s 9 and 10 would 
yield: 

k = f 3 k* + k (14) 

The results are shown in Figure 4, where the 
effective thermal conductivity is plotted versus sample 
average temperature for the three insulation densities. 
The experimentally measured apparent thermal 
conductivities are also presented in the figure. The 
error bars associated with the experimental data 
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presented in all the figures corresponded to 7.5% of 
reading uncertainty in accordance with the experimental 
uncertainty data presented in Table 1. The apparent 
thermal conductivity varies non-linearly with sample 
average temperature, increasing rapidly with increasing 
temperature due to the nonlinear nature of radiation 
heat transfer. This effect is more evident with the lower 
density insulation. As the density increases the solid 
conduction contribution increases, but the radiation heat 
transfer decreases more rapidly, resulting in a net 
decrease in the apparent thermal conductivity. The 
standard deviation between measured and calculated 
data for sample densities of 24, 48 and 96 kg/m was 8, 
6.3 and 5.4%, respectively. The agreement between the 
numerical and experimental results improved with 
increasing insulation density. 

For the sample density of 24 kg/m 3 the numerical 
data over-predicted the experimental results at lower 
temperatures. The contribution of solid conduction to 
the effective thermal conductivity decreases with 
decreasing insulation density in accordance with Eq. 

14. For the lowest density insulation solid conduction 
was almost negligible, and the effective thermal 
conductivity simply consisted of radiation. Therefore, 
the larger discrepancy between experimental and 
numerical results for the lowest density insulation was 
simply due to the radiation model. Because the values 
of optical thickness are low for this low-density 
insulation at low temperatures, the optically thick 
model begins to break down. However, at higher 
temperatures the optical thickness of the sample 
increased, resulting in closer agreement between 
experimental and numerical results. 

For sample densities of 48 and 96 kg/m 3 , the 
experimental and analytical results agree within the 
experimental accuracy. For sample density of 24 
kg/m 3 , agreement is not as good. If more accurate 
predictions are required an improved radiation model 
must be used. 

Thermal Conductivity in Presence of Gas 
Conduction 

The analytical model for gas conduction (Eq. 5a) 
and the two models for combined thermal conductivity 
( Eq. 9 and 10) were evaluated by comparison with 
experimental results with pressure values above 10 3 
torr. The value of the accomodation coefficient, a, in 
Eq. 5a was varied between zero and unity. The best 
results were obtained with a=l, and these are the 
results presented here. 

The variation of effective thermal conductivity 
with sample average temperature at environmental 
pressure of 100 torr for the three insulation samples is 
shown in Figure 5. The analytical predictions using 


both Eq.'s 9 and 10 were within the uncertainty range 
of experimental data for the three samples. The 
numerical results slightly under-predicted the 
experimental results at higher temperatures for sample 
densities of 24 and 48 kg/m 3 . This could attributed to 
uncertainties in gas conduction model and the coupling 
models used for combining solid and gas conduction. 

The variation of effective thermal conductivity 
with environmental pressure at nominal sample average 
temperatures of 50, 220, and 500°C for the three 
insulation densities is shown in Figures 6a through 6c. 
The increase in apparent thermal conductivity with 
increasing pressure is due to gas conduction. Gas 
conduction is negligible below 0.01 torr, rapidly 
increases between 0. 1 and 10 torr, and then 
asymptotically approaches a constant between 10 and 
760 torr. The standard deviations between measured 
and calculated data were 7.8 and 7.6% for the parallel 
model (Eq. 9) and the parallel/series model (Eq. 10), 
respectively. There was good agreement between 
analytical and experimental results for insulation 
density of 96 kg/m 3 over all temperatures and pressures, 
except for the data at 50°C at intermediate pressure 
levels (0.1 to 1 torr). The analytical results slightly 
under-predicted experimental results in this range. The 
parallel model (Eq. 9) matched the experimental results 
better for this insulation density. For the insulation 
density of 48 kg/m 3 , there was good agreement between 
numerical and experimental results over all pressures at 
sample average temperatures of 50 and 220°C. The 
numerical results under-predicted experimental results 
at pressures above 5 torr at 500°C, as discussed 
previously. For the insulation density of 24 kg/m 3 , the 
numerical results under-predicted experimental results 
at pressures above 10 torr at 500°C, and over-predicted 
at pressures below 0.01 torr at 220°C, as discussed 
previously. The parallel/series model (Eq. 10) was 
slightly better than the parallel model for insulation 
densities of 24 and 48 kg/m 3 . 

It can be concluded that the solid and gas 
conduction models and the two models for combined 
thermal conductivity produced satisfactory results. The 
optically thick assumption provides a good approximate 
technique for handling radiation heat transfer in fibrous 
insulation, with the accuracy of results improving with 
increasing insulation density. More detailed radiation 
models are required for more accurate modeling of 
fibrous insulation samples with density of 24 kg/m 3 . 
Insulation tested in metallic TPS under the X-33 
program 1 had a nominal density of 48 kg/m' and was 
several times as thick as the insulation samples tested in 
this study, resulting in a much higher optical thickness. 
Therefore, the optically thick assumption can be used 
for the analysis of its thermal performance. 
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Concluding Remarks 

The combined radiation and conduction heat 
transfer through fibrous insulation was solved 
numerically. The radiation heat transfer was modeled 
using the optically thick approximation. Various 
models for gas and solid conduction and their 
combination were evaluated. 

An experimental apparatus was designed and 
fabricated to measure the apparent thermal conductivity 
of high temperature insulation subject to large 
temperature differences representative of typical launch 
vehicle re-entry conditions. The insulation sample cold 
side was maintained near room temperature, while the 
hot side was heated to temperatures as high as 1000°C. 
The environmental pressure was varied from 1 x 10 4 to 
760 torr. The overall uncertainty for the apparent 
thermal conductivity measurements was between 5.5 to 
9.9%. The apparent thermal conductivity of Saffil™, 
an alumina fiber insulation, at nominal densities of 24, 
48, and 96 kg/m 3 was measured. 

The measured and predicted thermal 
conductivities were compared. It was found that the 
use of the optically thick approximation for modeling 
the radiation heat transfer through the insulation 
provided agreement within experimental accuracy over 
most temperatures and pressures. The agreement was 
not as good for the lowest density insulation sample 
studied. The models for gas and solid conduction and 
their coupling produced agreement within experimental 
accuracy. The best analytical model resulted in 
effective thermal conductivity predictions that were 
within 8% of experimental results for the tested range 
of pressures and temperatures. 
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Figure 1. Schematic diagram of thermal conductivity 
apparatus. 



Figure 2. Photograph of thermal conductivity apparatus 
with picture frame set on the water-cooled plate. 



Fig. 3a 



9 



Fig. 3b 

Figure 3. Contribution of different heat transfer 
mechanisms to local thermal conductivity for insulation 
density of 48 kg/m 1 with: a) P =10 4 torr, b) P=100 torr. 



Figure 4. Comparison of numerical and experimental 
thermal conductivity at 10 4 torr. 
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Figure 6. Comparison of numerical and experimental 
thermal conductivity versus pressure for insulation 
densities of a) 96, b) 48, c) 24 kg/m 3 . 


Figure 5. Comparison of numerical and experimental 
thermal conductivity at 100 torr. 
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